Non-linear spin to charge conversion in mesoscopic structures 
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Motivated by recent experiments [Vera-Marun et al., arXiv:1109.5969 , we formulate a non-linear 
theory of spin transport in quantum coherent conductors. We show how a mesoscopic constriction 
with energy-dependent transmission can convert a spin current injected by a spin accumulation into 
an electric signal, relying neither on magnetic nor exchange fields. When the transmission through 
the constriction is spin-independent, the spin-charge coupling is non-linear, with an electric signal 
that is quadratic in the accumulation. We estimate that gated mesoscopic constrictions have a 
sensitivity that allows to detect accumulations much smaller than a percent of the Fermi energy. 

PACS numbers: 72.25.Dc, 73.50.Fq, 75.76.+J, 73.23.-b 



Introduction. — Spin current detection and measure- 
ment protocols are for the most part based on ferro- 
magnetic contacts [H-Q ot Zeeman fields While 
efficient and well controlled, these schemes are not op- 
timized for miniaturization, because exchange and mag- 
netic fields have low spatial resolution and because they 
cannot detect the weak spin accumulations achievable 
in two dimensional electron gases such as GaAs het- 
erostructures, the platform of choice for sub-micron spin- 
tronics. To unleash the full potential of spintronics at 
the nanoscale, it is therefore imperative to find novel, 
all-electric protocols. In sub-micron structures, how- 
ever, reciprocity and other symmetry relations constrain 
the detection of spin currents in the linear response 
regime [H - 14 1 . In the absence of time-reversal symmetry 
breaking field and focusing on two-terminal geometries, 
these constraining rules can only be waived by going be- 
yond the linear response regime. At the nanoscale, this 
presents important theoretical challenges as local electric 
potentials must be determined self-consistently to ensure 
gauge invariance fl5j . 

In this manuscript we construct a mean-field non-linear 
theory of spin transport through sub-micron scale struc- 
tures. We use it to propose a protocol which converts 
the spin current injected by a spin accumulation into a 
charge signal via the energy-dependent transmission of 
a mesoscopic structure. While our scheme is general, 
we focus our discussion on quantum point contacts and 
Coulomb blockaded quantum dots, whose transmission 
is easily tunable by electric gate potentials. We show 
that the electric response is quadratic in the spin accu- 
mulation 5[i when the transmission is energy-dependent. 
A linear response arises only if the transmission is spin- 
dependent, which usually requires an external magnetic 
field. We foresee that our scheme has sufficient sensi- 
tivity to detect weak spin accumulations such as those 
that can be generated magnetoelectrically in GaAs het- 
erostructures jl6j], to which the magnetic spin detection 




FIG. 1: Measurement scheme: A mesoscopic constriction 
separates two non-magnetic terminals with spin accumu- 
lations <5/isi,2. The transmission coefficient through the 
constriction is gate tunable (voltage V g ), and when it is 
energy-dependent, an electric signal that is non- linear in 
SjJ, s i,2 arises. 

schemes are notoriously difficult to apply. 

We consider the standard measurement setup depicted 
in Fig. 1. One aims to detect the non-equilibrium spin 
accumulation drop below the two terminals, S^i s i ^ 5/i S 2, 
as an electric signal. The spin accumulation origin is not 
specified, be it ferromagnetic, magnetoelectric, or optical. 



A recent pioneering work, Ref. [17|, demonstrated that 
non-linear effects make the detection possible in graphene 
even without using ferromagnetic terminals, on which 
the linear Johnson-Silsbee method relies A volt- 

age quadratic in the spin accumulation arises due to the 
cncrgy-dcpendcncc of the graphene conductivity near the 
Dirac point. The bottom- line of our theory is that, while 
the drift-diffusion approach of Ref. [lTj is appropriate 
for bulk systems, it cannot be directly exported to sub- 
micron structures, where gauge invariance requires spe- 
cial care [l5j . Furthemore, unlike in graphene, the den- 
sity of states in GaAs heterostructures is mostly energy- 
independent, thus non- linear effects emerge only if fur- 
ther constrictions induce energy- dependent transmission 
T(E). The constriction, such as a Coulomb blockaded 
quantum dot, a resonant tunneling barrier, or a quan- 
tum point contact (QPC), is the active element in our 



2 



scheme, converting the spin to charge in proportion to 
d E T(E). This quantity, and this is the crucial point, is 
fully tunable electrically and independently of the spin 
accumulation itself, providing our method with versatil- 
ity necessary for practical spintronics. 

Theory calculation. -We model the detection circuit in 
Fig. 1 as a quantum scatterer connected to two elec- 
tron reservoirs, each with its own electrochemical poten- 
tial and spin accumulation, via two leads. We start by 
writing the current in lead i = 1, 2 in the spin subband 
a =t, I (alternatively a = ±1) [U] 

T i = lj dE^^S^S^-T^'iEMmffiE). 

w 

The transmission is the probability that a particle 
with spin a' injected from reservoir j exits the system 
with spin a into reservoir i. It depends on the particle 
energy E and the local electrostatic potential U(r). We 
consider a spin conserving transmission T aa oc S aa i with 



T u( E ) = T ( E ) + oST(E) 



(2) 



For a spin insensitive structure, the transmission differ- 
ence is zero, ST = 0, and T = T n = T u . Each lead is 
characterized by the number N? of transmission chan- 
nels it carries (whose weak energy- dependence we ne- 
glect), and the particle distribution f[(E) = f(E — /if), 
at the corresponding terminal, with the Fermi function 
f(x) = {cxp[(x — HF)/ksT] + 1} _1 . The electrochemi- 
cal potential of the spin subband a, measured from the 
Fermi energy \x E , is (e is the electron charge) 



/if = eVi + a <5/t s 



(3) 



where Vi is the applied voltage. 

Equation (JTJ) is current conserving due to the unitarity 
condition £V a = • To guarantee gauge invari- 
ance (i.e. that currents are invariant under an overall 
voltage shift) one has to take into account that U(r) is a 
function of the ap plie d voltages. Up to second order in 
/i's, the current is 15 
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Though the first linear term is explicitly gauge invariant, 
self-consistent conditions have to be imposed on U(r) in 
order to ensure that the non-linear terms also are gauge- 
invariant. The deviation SU (r) of the electrostatic profile 
from its equilibrium value, U(r) = U c(1 (r) + SU (r) , results 
from applied voltages and spin accumulations and the 
way they are injected into the scattering region. This 



deviation can be parametrized by characteristic poten- 
tials 15j, which generally speaking arc determined by 
self-consistent solutions to the Schrodingcr and Poisson 
equations. To restrain from these (here unnecessary) 
complications we neglect the spatial dependence of the 
potential changes, SU(r) = SU, and calculate the func- 
tional derivative of T[ 2 a using the identity 



J dr [8 u(r) T{ 



(E)} = -ed E T™ (E) 



(5) 



To solve for SU, we assume a symmetric probe, with 
equal, spin-independent coupling to both leads 



eSU = (/ij + ft 



t 



4)/4. 



(6) 



Using Eqs. (J2HH) wc get our main result, that the electri- 
cal current / 



l\ + 1\ is 



I = d e SV + G 2 {Sn 2 sl - + G 3 (<W - <W) 
+ G 4 (<fyt s i + Sns2) e SV. 



(7) 



The formula is explicitly gauge invariant, as it depends 
only on SV = V1—V2. The calculation is furthermore cur- 
rent conserving, with I2 = —I\ obtained by substituting 
Sfj, s i Sfj, S 2 and SV —> —SV. Wc sec the emergence of a 
non-linear spin to charge coupling term G-2(Sijl 2 s1 — S^i 2 ^), 
even in the absence of any rectification term oc SV 2 . That 
such a term is absent follows from our choice of a sym- 
metric potential SU, in agreement with Ref. 15j. There 



are four contributions to the current, with conductances 
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d E f)[d E ST(E)}, 


(8d) 



which we discuss in detail below, first for a spin- 
insensitive, second for a spin-sensitive constriction. 

Spin insensitive constriction. -"We first consider ST = 
in Eq. ((2J, in which case G3 = = G4, and focus 
our discussion on a gate-defined QPC in a 2DEG GaAs 
heterostructure, with energy-dependent transmission fl8j | 



T(E) = {l + cx P [-2tt(E - eaV g )/huj}Y 



(9) 



The transmission is easily tuned by an external gate 
voltage V g , with a sensitivity set by the QPC charac- 
teristic energy scale fvu) and a the "lever arm" convert- 
ing gate voltage into energy. We take typical values 
hit) = 180 /ieV, corresponding to the Zeeman energy of 
8 Tesla field at the g-factor g = —0.39, and a = 0.05. 
We further fix fip = 8 mcV, T = 0.1 K, and spin accu- 
mulations S[i S 2 = 0, and Sfisi = S/i s = 0.1% /xj?, which 
should be magnetoelectrically achievable flif . 
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FIG. 2: a) Non- linear conductance Gi given in Eq. (|8b[l . b- 
d) Current, Eq. © for zero bias, 8V=0, as a function of b) 
gate voltage determining the QPC transmission [See Eq. ©], 
c) temperature and d) QPC energy resolution in units of the 
magnetic field. 

We are now ready to investigate the electric response of 
the circuit. First, we assume that both leads are held at 
the same potential. Even in this case, the spin accumula- 
tion generates a finite current, due to the second term in 
Eq. (|7|). Its magnitude is determined by the non-linear 
conductance Gi in Eq. (|8b[) . which we plot in Fig. 2a. 
It is proportional to 5eT, and thus maximal when the 
QPC is half open, T = 0.5. The current is then 

I = G ^A %\ /9 M , (io) 

and we plot it in Fig. 2b. For our choice of parame- 
ters, the current is of the order of tens of pA, which 
is well above the experimental detection limit. The de- 
pendence of the current signal on max(fesT, Hljj/2tt) is 
demonstrated in Fig. 2c and d. Decreasing ksT at fixed 
Huj (Hu at fixed fceT), the signal first increases before it 
saturates when fcflT~ fiw. 

Alternatively, terminal 2 can be operated as a floating 
probe. In this case, a finite voltage drop develops, which 
we find by setting I\ = 0. One obtains 

eSV = -G*6ii 2 JG x . (11) 

We see that the current is converted into a voltage by the 
linear conductance G\, given in Eq. (|8a|) and plotted in 
Fig. 3a. This agrees with the already mentioned absence 
of a rectification term in our symmetric QPC We 
plot the signal voltage in Fig. 3b. In the region where the 
QPC is closed (V g -C — Huj), Eq. (|TT|) gives an unphysical 
saturation of SV (dashed line). To remove this artifact, 
we add a small constant to Gi, which enforces that Sfi s i 
does not influence V2 if the QPC is closed. Then the 




FIG. 3: a) Linear conductance Gi, Eq. (|8a[) . b) Voltage, 
calculated according to Eq. (|lll) (dashed line) and adding a 
constant 0.1 e 2 /h to Gi (solid line), c) Current in external 
magnetic field antiparallel (solid line) and parallel (dashed 
line) to the spin polarization direction. The arrow denotes 
the zero current position, the thin line is a guide to the eye. 
d) Second derivative of the current with respect to V and B. 

electric signals (current or voltage), in the two protocols 
just discussed behave similarly, I,V oc G2<5/U 2 . For spin- 
insensitive constrictions, we see that the electric response 
is quadratic in the spin accumulation. 

This result is not specific to a QPC, which we next 
replace by a Coulomb-blockaded quantum dot. Ne- 
glecting inelastic processes and near resonance, its low- 
temperature transmission is given by [19| T(E) = 
T^Y^ /[{E-E Q f /h 2 +{T /2)% with the tunneling rates 
rW and r( 2 ) of the resonant level to the left and right 
leads, r = rW + T^ 2 ' and the resonance peak position 
E . At E - E Q = ±hT/2y/Z, d E T{E) takes its maximal 
value iOrWr^/^r 3 . Because T(E) differs from its 
QPC expression, Eq. ©, in its energy-dependence, the 
shape of the quantities plotted in Fig. [5] will be different; 
most notably, the signal changes sign upon crossing the 
resonance. However, the maximal current magnitude is 
still given by Eq. (jlOp with T replacing u>. 

Spin sensitive constriction— We next consider the case 
when ST 7^ in Eq. ©, when the QPC is made spin 
sensitive, e.g. by an external Zeeman field. We assume 
that the field is parallel to the spin accumulation S[m s i 
and that it is sufficiently weak that the latter is not in- 
fluenced. Equation (j9]) becomes 

T aa (E) = {1 + cxp[-2n(E - vpB - eaVg)/^]}' 1 , 

(12) 

where the Zeeman energy [iB is added/subtracted from 
the electron's energy depending on its spin. Here, B is 
the magnetic field, /i = (g/2)fis and (is is the Bohr 
magneton. To linear order in B, we then have ST = 



4 



—fiBdET. We see that a term linear in the spin accu- 
mulation has appeared, whose magnitude is given by the 
conductance G3, Eq. ([5c]) . This is a linear response term, 
similar to the one reported in Ref. giving an odd 



magnetoresponse of the electric current through a QPC 
in the presence of a spin current. This term allows to find 
the sign of the spin accumulation, which is impossible for 
a spin insensitive probe. This is demonstrated in Fig. 3c, 
where we plot the current as a function of the magnetic 
field. Applying the field antiparallel (for g < as in GaAs 
heterostructures) to the spin polarization direction, the 
Zeeman energy penalty compensates for larger transmis- 
sion at higher energies. The compensation is exact (the 
current becomes zero), if 



8(i s = gfiBB c 



(13) 



Remarkably, determining the compensation field B c 
alone allows to measure both the magnitude and direc- 
tion of the spin accumulation. 

For spin-sensitive constrictions, the conductance G4 
gives an interesting contribution to the current, which is 
coupled to the average spin accumulation 5/i s i+<fyi S 2 and 
the voltage bias 8V. We rewrite this contribution as 



Gi(Sfx 2 + Sfx 1 )eSV = G 4 



'Mi 



t 
M2 



4. 



eSV. 



(14) 

which makes it clear that this term describes two differ- 
ent rectification currents in the two spin subbands which 
are uncompensated if (i) the transport in the subbands 
happens at different energies, (ii) there is a finite bias 
and (iii) the probe transmission is both spin and energy 
sensitive. Experimentally, this contribution can be iden- 
tified from the current derivative with respect to both 
the applied bias and the external field, as in Ref. poj ). 
We plot this contribution in Fig. 3d where the antisym- 
metric shape of G4 oc d\T , strongly contrasts with the 
symmetric conductances G2, and G3. 

Conclusions. -We have shown how spin accumulations 
can be converted into electric signals in mesoscopic sys- 
tems with energy- dependent, but spin-conserving trans- 
mission. When transport is spin-independent and in the 
absence of voltage bias, the conversion occurs in the non- 
linear regime, and the electric signal is quadratic in the 
spin accumulation. In sub-micron structures, such non- 
linearities have to be treated self-consistently in local 
electrostatic potentials generated by the finite applied 
biases. We did that within a simplified mean-field ap- 
proach, which resides in neglecting the spatial structure 
of the potential changes SU(r). It is reassuring that the 
non-linear signal arises within this approximation, which 
restrains from details of the constriction. Further, device 
specific, spin rectifying effects may arise from the spa- 
tial effects in <5t/(r), along the lines of Refs. [IH We also 
note that for no applied bias our approximation becomes 
exact since 8U — s- 0. 



In the case of transmission through a QPC, Eq. (p~0|) 
suggests that it is its energy resolution, with hoj ~ 2 — 
3K typically, rather than the temperature, which limits 
the signal magnitude and hence the sensitivity of our 
approach. Alternatively, a Coulomb-blockaded quantum 
dot can be used, where the resolution is given by the 
tunneling width which can easily reach hT ~ 0.1 K or 
less (see e.g. (HJ). Close to resonance, we would expect 
such a quantum dot to enhance the signal by at least 
one order of magnitude compared to the results shown in 
Fig. 2d. 
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